Entanglement of finite cyclic chains at factorizing fields 
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We examine the entanglement of cyclic spin 1/2 chains with anisotropic XYZ Heisenberg cou- 
plings of arbitrary range at transverse factorizing magnetic fields. At these fields the system exhibits 
a degenerate symmetry-breaking separable ground state (GS). It is shown, however, that the side 
limits of the GS pairwise entanglement at these fields are actually non-zero in finite chains, corre- 
sponding such fields to a GS spin-parity transition. These limits exhibit universal properties like 
being independent of the pair separation and interaction range, and are directly related to the mag- 
netization jump. Illustrative exact results are shown for chains with I) full range and II) nearest 
neighbor couplings. Global entanglement properties at such points are also discussed. 

PACS numbers: 03.67.Mn, 03.65.Ud, 75.10.Jm 



Quantum entanglement is well recognized as a funda- 
mental resource in quantum information science It 
also provides a new perspective for the analysis of quan- 
tum many-body systems, allowing to identify the gen- 
uine quantum correlations [2-4]. An important result 
for quantum spin chains with finite range interactions 
is that in contrast with the correlation length, the pair- 
wise entanglement range does not necessarily diverge at 
a quantum phase transition For instance, it remains 
confined to just first and second neighbors in a near- 
est neighbor Ising chain placed in a transverse magnetic 
field [2j. It can, however, diverge at a different point. 
Spin chains with anisotropic coupling exhibit a remark- 
able factorizing field 0-0] , where the system possesses a 
separable GS and hence entanglement vanishes in princi- 
ple, although it was shown to reach infinite range in its 
vicinity Q- 

Previous analyses were focused on large systems. The 
purpose of this work is to investigate the entanglement of 
finite chains, relevant for quantum information process- 
ing, with anisotropic XYZ coupling of arbitrary range, 
exactly at transverse factorizing fields. Although it may 
seem that any type of entanglement will vanish at such 
points, it should be noticed that separable ground states 
break a fundamental symmetry of these chains (the S z 
parity or global phase flip [2J) and are hence degenerate 
(and non-orthogonal), the factorizing field correspond- 
ing actually to a GS transition between opposite parity 
states. As a consequence, pairwise entanglement in finite 
chains will be shown to approach distinct non-zero side 
limits at the factorizing field, which do not depend on the 
pair separation or coupling range and whose average is 
directly measurable through the magnitude of the magne- 
tization jump at the transition. Moreover, even the pro- 
jector onto the GS subspace remains entangled at these 
fields. Although these effects become negligible in large 
anisotropic chains, where opposite parity ground states 
are nearly degenerate, they will be shown to be quite 
prominent in small finite chains and, moreover, to remain 
appreciable for increasing sizes if the XY anisotropy be- 
comes sufficiently small. Present results are therefore 
particularly relevant for chains close to the XXZ limit. 
We derive first general exact results valid for any range, 



describing then illustrative exact results for chains with 
full range and nearest neighbor couplings. The exact 
results in the last case are obtained through the Jordan- 
Wigner mapping and its analytic parity dependent diag- 
onalization. 

We consider a cyclic chain of n qubits or spins interact- 
ing through an XYZ Heisenberg coupling with arbitrary 
common range in a transverse magnetic field b. Denoting 
with s l the spin at site i, the Hamiltonian reads 



H = bS z -J2 r j-i( v xSls J x + Vys^l + v z s\sl) 



(1) 



i<3 



= bS z -^2 r j-i[U v + s + s - + v - sl + s + + h.c) + v z s l z s J z ] , 

where S = " =1 s\ v ± = ( v x ^ v v)l^ an d r l = r n-i 
for I = j — i = 1, . . . , n — 1. Without loss of generality 
we can here assume b > and v x > \v y \ (i.e., v± > 0), 
with r\ arbitrary. We will be interested in the attractive 
(ferromagnetic) case 77 > V I, although the following 
considerations are general. Since H conserves the S z - 
parity, 



[H, P z ] =0, P z = exp[i7r(S 2 + n/2)] , 



(2) 



its nondegcncratc cigcnstates will have definite parity 

Pz = ±1. 

Let us now examine the conditions for which a com- 
pletely symmetric separable state of the form 

n 

\9) = [](cos \9\U) + sin i0|t;)) = exp[i6S y }\0) , (3) 

i=l 

where s*|J^) = — h\Jri) an d |0) = IlilW) can be an 
exact eigenstate of (Q}. This state is fully aligned along 
an axes z' forming an angle 6 with the z axes, such that 
S Z '\0) = —hn\9), breaking parity symmetry for £ (0,7r). 
Replacing s\ x = s\, x , cos8±s l x , z , sin# in (TTJ), it is easily 
seen that these conditions are 



cos(9 = ±^x, x 



_ u y 



b = r(v x — v z ) cost 



n-l 

1=1 



(4) 
(5) 



2 



where Eq. ([5]) is required for 9 € (0,tt), i.e., x G [0, 1) (in 
the XX Z case x = 1 (v y = v x ) both |0) and \tt) are trivial 
eigenstates of ([T]) for all fields &). Such parity breaking 
separable eigenstate is then feasible for x <= [0,1) (i.e., 



< v y < v x if v x > \v y \) and b = ±6 S , with 
b s ee - v z )y/x, 



(6) 



the factorizing field. The state ([3]) will depend on the 
anisotropy x but not on the factors ri, being then inde- 
pendent of the interaction range. 

It is also apparent that both \9) and |— 6) = P z \9) are 
degenerate eigenstates of H at b = b Sl with energy 



(9\H\9) 



-^n[bcos9 i 
■\nr{v x 



v 



\r{v x sin 2 I 
v z ) . 



+ V z COS' 



0)] 



(7) 



Hence, at b — ±6 S two levels of opposite parity necessar- 
ily cross, enabling the formation of these eigenstates. 

Let us remark that in the attractive case r; > VZ, the 
state that minimizes (H) among separable states (i.e., the 
mean field approximate GS) is precisely of the form ([3]) 
V&, with \9\ determined by Eq. if |&| < b c = r(v x — v z ) 
(parity-breaking solution) and 9 — otherwise. Hence, 
in this case the factorizing field can be seen as that where 
the mean field GS becomes an exact eigenstate (i.e., the 
exact GS, as shown below). 

The states |± 9) will then form a basis of the corre- 
sponding eigenspace at b — b s (assumed of dimension 2), 
which is non- orthogonal for 9 ^ ir/2: (—9\9) = cos™ 9. A 
proper orthonormal basis conserving parity symmetry is 
provided by the entangled states 



\0±) = 



±1- 



y/2(l ± COS 

E 



V2sm k % cos"- fc | ., , 
- - A S^|0), 



kly/1 ± cos" 9 



(8a) 
(8b) 



which satisfy P z \9±) — ±|0±) and are the actual eigen- 
states of H in each parity subspace at b — b s . These 
states (and not the states | ± 9)) are the actual limits of 
the corresponding exact eigenstates | X P ± (6)) (which have 
definite parity) for b — > b s . 

In the attractive case n > V I (with \v y \ < v x ), the 
states \9±) (and hence |± 9)) are ground states of H at 
b = b s : The exact GS ^^(b)) in each parity subspace 
must have expansion coefficients all of the same sign in 
the standard computational basis (i.e., that of separable 
states with definite values of {s z }) in order to minimize 
the average energy, since the average of the off-diagonal 
XY term in (fT]) can only increase (or eventually stay 
constant) for different signs (as rj_i > 0, v± > 0) while 
those of the diagonal terms bS z and v z s l z s 3 z are sign inde- 
pendent. Hence, |*q (6)) cannot be orthogonal to \9±), 
whose expansion coefficients in this basis are all non-zero 
and of the same sign (Eq. (|8b[) ). and must then coincide 
with \9±) at b = b s . 

Thus, in the attractive case \9±) represent the side lim- 
its lim b ^ b ± \^a(b)) of the exact GS \^n{b)) in the whole 



space at b = b s , which undergoes there a \9-) —> \9 + ) 
parity transition (actually the last parity transition as b 
increases, as will be shown in the examples). 

The pairwise entanglement in the states \9±) depends 
essentially on the overlap (—9\9). When orthogonal 
(9 = 7r/2), they are generalized GHZ states which, al- 
though globally entangled, exhibit no pairwise entangle- 
ment (for n > 2). Moreover, in this case the normalized 
projector onto the space spanned by the states \9±), 



po = h(\e+)(e+\ + \e-){0-\) 



(9) 



which represents in the attractive case the T — > limit of 
the thermal mixed state p(T) cx exp[— H/kT] at b — b s , is 
fully separable (i.e., a convex combination of projectors 
onto separable states) as po = h(\0){6\+\-6){— 9\) for 9 = 
it/2. In contrast, for 9 E (0,7r/2) both states ([8]) as well 
as the mixed state ^ will be shown to exhibit a uniform 
non-zero entanglement between any two spins (note that 
the projector onto this subspace is no longer the sum of 
the individual projectors | ± 9}{±9\ when (— 9\9) ^ 0). 

Let us first evaluate the pairwise concurrence @ (a 
measure of pairwise entanglement) in the states \9±). As 
a consequence of ^ and the cyclic nature of H, the re- 
duced two spin density matrix pij in any non-degenerate 
eigenstate or in p(T) , will commute with the reduced par- 
ity P« = e "r(sl+sl+i) and depend just on I = \i-j\. The 
ensuing concurrence C; ee C(pij) takes the form 



Ci = 2Max[|( 



-01,0] 



(10) 



where af = (s 4 + s 3 ± ), pi = \- (s l z sl), qi = [(± - Pl ) 2 - 
(si) 2 } 1 / 2 . If \a+\ > Pl [\a-\ > qi ) C, is of even (odd) 
parity type, i.e., parallel (antiparallel) Q, as in Bell state 
« ITT) + (It-i-) + lit))- Just one of these inequalities 
can be satisfied in a given state. 

In the states ([8]), a\ — -jsin 2 ^^, pi = a^, (s l z ) = 

-\ cos6>7^, with 7^ = l± Z™f s ~l 6 and v = ±. We then 
obtain Q{\9±)) = C± V I, with (assuming 9 € (0,tt/2]) 

cos™" 2 9 



C± = sin' 



1 ± cos™ I 



(i-x) 



X 



n/2-1 



1 ± x™ /2 



(11a) 
(lib) 



Thus, C_ > C+ > 0, with C+ (C_) parallel (antiparallel). 
Note that for 9 ->• (x 1), 

C-^2/n, 

as in this limit \9+) |0) but \9-) -4- |1) = -^S+ 10), 
which is an W -state Q (2/n is in fact the maximum value 
that can be attained by the concurrence in fully symmet- 
ric states [HI). As 9 increases, C_ decreases while C+ 
becomes maximum at 9 ~ 1.6 /y/n (see Eq.ll3p. vanishing 
both for 9 tt/2 (x 0) if n > 2. 

In the attractive case the values (jllb[) represent the 
universal side limits C± = lim^^i Ci(b) of the GS con- 
currences Ci(b) at b = b s , valid for any separation I or 



3 



For x — > 1 they correctly approach 
|1) — > [0) transition taking place at b = b c 
limit [ill 



interaction range, 
those for the 

in the XX Z limit [11] (where b s -t b c ). 

The concurrence jump C-— C+ determines, noticeably, 
the concurrence Co = Ci(po) in the GS mixture ([9j, 



Co = §(CL - C + ) = (1 - X )j 



71-1 



X" 



(12) 



(see also Eq. IT71) . which is of antiparallel type. It is a 
decreasing function of 0, starting at 1/n for 6—^-0. In the 
attractive case, Eq. (fT2")l represents the common T — s- 
limit of the thermal concurrences C; (T) at b — b s for any 
separation I and coupling range. 

Although for fixed x < 1, C± become exponentially 
small as n increases, the rescaled concurrences nC± re- 
main finite for small anisotropy % = 1 — S/n. For large n 
and fixed <5, we obtain from (fTTj) -(fT2 j) the n-independent 
limits 



c± = nC± 
c = nC ' 



8e- 5 / 2 /{l±e- s ' 2 ), 



6e-'/(l - e~% 



-6\ 



(13) 
(14) 



depicted in Fig. Q] While c_ and c = c_(2<5)/2 arc 
decreasing functions of 5, c + is maximum at 5 = 2[1 + 
u>(e -1 )] ss 2.56, where c+ = 2u>(e~ 1 ) s» 0.56 (w(x) is the 
productlog function, such that x — we w ). We note also 
that Co > c+ for S < 2 In 2. 

The mean rescaled concurrence (c++c_)/2 determines, 
remarkably, the total magnetization jump at b = b s : 



AM = (6-\S z \6-) - (8+\S z \6 + ) 



nsm 2 0- 



cos 



1 



cos^ 



i(c++c_)^, (15a) 



which represents as well the slope of the energy gap AE w 
(b — b s )AM between the odd and even GS at b — b s . For 
large n and fixed 6, 



AM « (c+ + c_)/2 = 5e- s/2 /(l - er s ) , 



(15b) 



remaining finite and providing a direct way to determine 
the average rescaled concurrence at b s . 

As illustration, we now show exact results for the 
concurrence in I) a fully connected chain with constant 
ri 12] and II) a chain with nearest neighbor coupling 
(n = <5;,i + <$;,n-i)- In I) we set n = 2/(n— 1) V I such 
that r = 1 in I and II (Eq. [5]) . The factorizing field © 
and the energy ([7]) are then the same in I and II for fixed 
v x,y,z- We will consider v x > and v z — (XY case). 

In I, the GS can be obtained numerically by diagonal- 
izing H in the subspace of maximum total spin states (to 
which it belongs) as [H, S 2 } = 0: 



Hj 



bS z 



E v,{Sl-\n)/{n-l). 



H=x,y 



The fixed parity GS is then of the form X)fc ovon w kS^_\0), 
leading to I independent elements = (S+)/c n , — 




FIG. 1. (Color online) Ground state rescaled concurrences 
at the factorizing field vs. scaled anisotropy parameter 8. c± 
denote the side limits (|13[1 . cq the value at b — b a (Eq. 1 14)1 . 
color indicating the antiparallel (c_, Co) or parallel (c+) type. 
The dotted line depicts the magnetization jump (|15[1 . These 
curves hold for any spin pair and Hamiltonian of the form {!]). 



(n 2 /4- (S 2 ))/c„, {s\s\+ l ) = {{Si) - n/4)/c„, with c„ = 
n(n — 1). 

In II, the Hamiltonian can be solved analytically for 
any finite n by means of the Jordan- Wigner transforma- 
tion [l3|], which allows to rewrite H, for each value (±) 
of the parity P z , as a quadratic form in fermion operators 
c\, a defined by cj = sf exp[-i7r s j~ s 7l : 

n 
i=l 

= E A fc(4«fc-s)» (16) 

fee-RT ± 

where n + 1 = 1, T]^ = 1, rj^ = 1 — 2dm and 

A 2 = {b — u + coswfc) 2 + v 2 _ sin 2 w^, = 2nk/n , 

with A+ = A_ = {0,...,n-l} 

(i.e., k half-integer (integer) for positive (negative) par- 
ity). The diagonal form (fir?)) is obtained through a 
discrete parity- dependent Fourier transform 11] cj 



^fcei<"± e "" fcJ ' c k) followed by a BCS transforma- 
tion c'J = u k a\ + v k a n ^ k , c' n _ k = u k a n ^ k - v k a\ to 
quasiparticle fermion operators a k , at, with u 2 ,i> 2 = 
|[1 ± (b — i;+ cos Wfe)/Afc]. For b > we set Afc > 
for fc 7^ and Ao = v + — 6, such that the quasiparticle 
vacuum in HJ I is odd and the lowest energies for each 
parity are Efj = - \ J2 k £ k± At b = b s = ^Jn^ (Eq. 

©), X k = v+ — 6 s coswfe and Efj = E7j = —nv+/2, in 
full agreement with Eq. ([7]). 

The concurrences in the fixed parity GS can be ob- 
tained from the contractions /; = {c\cj)± — \o~ij, gi = 
(c\cj)± and the use of Wick's theorem [lj|, leading 

to (4) = fo, (44) - fS - ft + 9? and a t = 

j[det{A+) Tdet{A t )], with (Af)ij = 2(/i_ i±1 +gi-j±i) 
I x I matrices. 

In both I and II, as b increases from + , [n/2] GS 
parity transitions ± — > =F take place if X € (0: 1] ( as in 
the XX Z case (ll|), the last one (— —> +) at b = b s . 
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FIG. 2. (Color online) Ground state rescaled concurrences c; 
vs. magnetic field b for n = 10 qubits and anisotropy x = 
1 — S/n — 0.75 (left) and 0.5 (right). Top panels correspond 
to a fully connected chain (case I), where ci is the same for 
all separations I, bottom panels to nearest neighbor coupling 
(case II), where ci is shown for all separations. Horizontal 
dotted lines indicate the limit values (|11[) and (|12[) . At b = b s , 
ci changes from antiparallel (red) to parallel (blue) type, the 
side limits being non-zero and identical V I, and the same in 
I and II. Large dots indicate concurrence values at the n/2 
parity transitions (for the mixture of both ground states), 
with that at b — b s given by Eq. (fl2|) V I in I and II. 



They are clearly visible for low n if x is not small, i.e., 
if 6 = n(l — x) is n °t too large, as seen in Fig. [2] for 
n = 10 qubits. For b — > b^, all GS concurrences Of 1 are 
seen to approach the same side limits (fTTj) in both I and 
II, which are non- negligible. For 6 = 2.5, Or reaches in 
fact its maximum for b — > b s ~ 0.87v x in I, and also in II 
if Z > 2. For 5 — 5 the side limits are still noticeable but 
nearly coincident, implying a negligible cq (Eg. I12[) . 

The behavior of cj~ for n — 50 qubits at the same values 
of S (now x — 0.95 and 0.9), depicted in Fig. El is seen 
to be the same as in Fig. [5] in the vicinity of b s . All pairs 
become entangled as b — > b s , with Ci approaching the 
common values ([lip in I and II, now well approximated 
by Eqs. (HID. 

Within each parity subspace, the factorizing field is dis- 
tinguished as that where all GS concurrences Cf~ cross 
at the values (fTTj) (rather than vanish) , as seen in the top 
panel of Fig. 21 Moreover, in II the ordering of concur- 
rences C, becomes inverted at b - 



0.96 



3s 0.98 b/vx i.o 



0.94 



bs 0.96 0.98 b/Vx 



b s : C ; ~ (C ; + ) increases 
with increasing separation I for b just above (below) b Sl 
as C; (6) is linear close to b s . Note that Cf(b) vanishes 
at a lower field b^ < b s , becoming antiparallel for b < b^ . 

On the other hand, at sufficiently low temperatures, b s 
can be identified as the field where all thermal concur- 
rences Ci(T) cross at the value (fT2"]l . as seen in the bot- 
tom panel of Fig. @] C;(T) vanishes at a slightly larger 
Z-dependent field bi(T) > b s , remaining antiparallel un- 



FIG. 3. (Color online) Same details as Fig. [2] for n = 50 
qubits. Bottom panels depict results in the vicinity of b s , the 
inset those in the full interval. The limit values at b = b s are 
again the same for all separations and identical in I and II. 



til bi(T). To understand this effect, we note that in a 
general mixture 

p g = q\e+){0 + \ + (l-q)\0-)(e-\, qe [0,1] 

the concurrence C(q) = Ci{p q ) is 



C(q) = \l-q/q c \C-, q e = ±(l + < 



(17) 



which generalizes Eq. (TT21 (recovered for q = 1/2). C(q) 
is antiparallel (parallel) for q < q c (> q c ) and zero at 
q = q c > 1/2, where p q becomes completely separa- 
ble (p qe = |(|0}(0| + |- 0){-6\)). Separability requires 
then a slightly greater weight in \9 + ) due to its lower 
concurrence. Hence, at low T > Ci(T) vanishes and 
changes from antiparallel to parallel at a slightly larger 
field bi(T) > b s , where the positive parity GS has a 
higher weight in the thermal mixture. In case II this 
entails the surprising result that in the narrow interval 
b s < b < bi(T), the thermal concurrence C;(T) will in- 
crease with increasing I, since it is still driven by the 
lowest odd state V I. 

Let us finally mention that for any chain partition 
(L,n — L), the Schmidt number of the states ([8]) is 2. 
Their Schmidt decomposition [l| is 



(18) 



Pl± = 




(l + l/COS L 0)(1±2/COS 

2(1 ± cos" 6) 



, (i/ = ±)(19) 



where |6>£) = (\6 L ) ± |- 6> L ))/ 1 /2(1 ± cos L 6) denotes the 
analogous fixed parity states for L spins and p u L± the 
eigenvalues of the ensuing reduced density. The entan- 
glement between L and n — L spins can be measured 
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0.972 b s b/v x 0-978 



FIG. 4. (Color online) Top: Rescaled concurrences c; in the 
lowest state for odd (upper curves, in red) and even (lower 
curves, in blue) parity P, for case II with 5 — 2.5 and n = 50. 
Numbers indicate the separation I. At b = b s an inver- 
sion in the ordering of the c;'s with / takes place. Dashed 
lines depict the antiparallel concurrence in the even states. 
Bottom: Thermal concurrences in the previous system at 
kT = 5 x lO -4 ^ (solid lines) and at T = (dotted lines). 

through the entropy Sf; = - Y, v Pl± 1o S 2 Pl± or equiva- 
lently, the "global" concurrence Cf — y/2(l — J2u(Pl) 2 ) 
[l4| (square root of the tangle for a pure state), which is 
just an increasing function of Sf (Sf,C^ G [0, 1]): 

u l~ z \/Pl±Pl±- 1±y™/2 

where we have replaced cos6> = y/x- Hence, > C~£, 
with C^; increasing functions of 9, i.e., decreasing func- 



tions of x, in contrast with the pairwise concurrences C±. 
For 9 -+ n/2 ( X -> 0), Cf -> 1 (GHZ limit), whereas for 
9 -> ( X -> 1), Cl -> but 



Cl -> 2y/L(n- L)/n, 

(W-state limit), in which case ~ (L/n)[l— log 2 (i/n)] 
for L -C n. Thus, within the bounds imposed by a 
Schmidt number 2, the behavior of and with i 
is "non-critical" (i.e. saturated) for low \ (large 5) and 
"critical" (non-saturated) for x 1 (l° w ^) and nega- 
tive parity. It is also verified that Cf~ > y/n — 1C± (in 
agreement with the general inequality C£ =1 > J2?=i @i 
saturation reached for 9 — > 0, where C^/C± » 
y^T[l + I^ 2 (n-2)]. 

In summary, we have shown that due to the 5 Z parity 
conservation, the GS of finite cyclic chains with attrac- 
tive couplings of the form ([1} remains entangled as the 
factorizing field b s is approached, undergoing at b s the 
last parity transition and exhibiting for b — > bf universal 
entanglement properties, "intermediate" between those 
of GHZ and W-states. This field plays thus the role of a 
"quantum critical field" for small chains, with the pair- 
wise concurrence reaching infinite range and approaching 
distinct side limits which are independent of the pair sep- 
aration and interaction range. Their average is directly 
measurable through the GS magnetization jump (|15al) . 
which provides then a signature of the present effects, 
while their difference determines the concurrence of the 
GS mixture (|12l) . These effects remain appreciable for 
increasing n if the anisotropy becomes sufficiently small 
(finite S), i.e., for chains close to the XXZ limit. More- 
over, within a fixed parity subspace (and also at suffi- 
ciently low T > 0), b s is singled out as the field where 
all pairwise concurrences cross, the ordering with separa- 
tion becoming inverted as b crosses b s . Type, range and 
even ordering of the pairwise entanglement can thus be 
controlled by tuning the field around b s . 

The authors acknowledge support from CIC (RR) and 
CONICET (NC,JMM) of Argentina. 
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